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Abstract 

We compute the proton structure function at small x and large Q 2 at next-to- 
leading order in a s (Q 2 ), including summations of all leading and subleading logarithms 
of Q 2 and 1/x in a way consistent with momentum conservation. We perform a detailed 
comparison to the 1993 HERA data, and show that they may be used to determine 
a s (M|) = 0.120 ± 0.005(exp) ± 0.009(th). The theoretical error is dominated by the 
renormalization and factorization scheme ambiguities. 
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Experimental knowledge of the proton structure function F<i at large Q 2 and small x 
has dramatically improved since data from HERA have become available Whereas 
the data are in beautiful agreement || with the leading-order prediction of perturbative 
QCD 0], and in particular display the double scaling behaviour which follows || from it, 
the improvement in experimental accuracy now opens up the possibility of more detailed 
tests of the theory. Indeed, reliable perturbative calculations in the small x region turn out 
to be possible thanks to the fact that perturbation theory can be reorganized to give an 
evolution equation which provides order by order a convergent sum of all leading logarithms 
of x and Q 2 ||. Here we will show that full next-to leading order (NLO) computations 
give an excellent description of all the HERA data. Furthermore, this description is largely 
independent of any detailed nonperturbative (thus uncalculable) information, such as the 
detailed shape of input parton distributions, but is relatively sensitive to the precise value 
of the strong coupling a s , which may thus be determined using presently available data 0. 

The determination of a s from the shape of F^ at small x which we will discuss 
here is similar to that used in the analysis of large x structure functions, however it 
also shares some of the features of measurements based on less inclusive quantities (see 
ref. In fact, many of the more inclusive measurements (for example e + e~ — > hadrons) 
require precise data with very high statistics, since it is necessary to measure the small 
<5 2 -dependent deviations from an a s -independent leading order result. The determinations 
from scaling violations in structure functions obtained from fixed target deep inelastic 
scattering are even more difficult, since there the leading behaviour (the parton distribution 
functions at some initial scale) cannot be reliably computed, but must be fitted. By 
contrast, measurements based on less inclusive quantities which are essentially proportional 
to a s at leading order (for example jet rates, event shapes and energy correlations) can 
be made with much lower statistics: the dominant error then comes from hadronization 
uncertainties. In this respect the determination of a s from at small x discussed here 
is akin to the jet determination, and, in particular, it can be made with relatively low 
statistics; however here there are no hadronization effects, and higher twist corrections 
will turn out to be relatively small. 

The possibility of determining a s from the presently available HERA data P],H for 
F% (x, Q 2 ) follows directly from the empirical fact that the data exhibit double scaling |J: 
throughout the HERA kinematic region In F^ rises linearly in a = y^C, but is independent 
of p = a/£/C (where, as in ref.0, £ = In £ = In ^, t = lnQ 2 /A 2 ). More precisely, 
RfF2, where Rf is proportional to cr 1 / 2 / oexp(— 27a + Sa/p), is independent of both a and 
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p. This may be seen immediately from the scaling plots fig. [l|. The two parameters 7 

At, 
0o 



a/12//3q, 5 = (11 + ~^f) I Po depend inversely on the rate of running of a s (t) = 4^ + 0(jz); 



(3q = 11 — I nf is the leading coefficient of the /3-function. The fact that the HERA data 
exhibit double scaling is thus direct evidence for the running of a s (t) ||. Furthermore, 
the growth of F2 at small x and large Q 2 is due primarily to the collinear singularity in 
the triple gluon vertex (at leading order in perturbation theory), and it is thus directly 
and simply related to the actual value of a s . This is seen most clearly by rewriting the 
predicted asymptotic leading order behaviour of Ff as 



Since the third term rises more slowly than any power of a s , at very large t (and fixed 
x) the second term dominates and Bjorken scaling is eventually recovered. However at 
smaller values of t and small enough values of x (and in particularly in the HERA region 
H) the strong growth in the third term overwhelms the power fall-off in the second, and 
F2 rises quickly as a s falls. It is this strong dependence of the leading behaviour of F2 at 
small x on the value of a s which makes it possible to make a precise determination of a s 
using the existing HERA data Indeed, fitting the functional form (^ to the data, 

fixing Qo = 1 GeV and xq = 0.1, as in 0, gives a s (Mz) = 0.125 ± 0.002 (at one loop), 
with a x 2 °f 68/121, which is at least encouraging. 

Of course, to actually determine a s requires at least a full NLO computation, in order 
to fix the scale. Indeed, a two loop computation |§ reveals that NLO effects lead to a 
reduction in rate of growth of F2 in the HERA region: although the asymptotic slope 27 
is approached at very large cr, subasymptotically the slope is somewhat lower. This is 
essentially due to the fact that the leading singularity in the two-loop splitting function 



Pf is of opposite sign to that at leading order: in MS scheme 

P«(x, t) = \2C A - ^ ~ 12C ^ SM] 1 + . . , (2) 

9 2n J x 

where Ca = 3, Cp = |, Tr = | for N c = 3. The slope of the rise in F 2 is accordingly 
reduced. Indeed this reduction can now be seen in the data (see table 4 and fig. |l|): 
agreement with the data is significantly improved, essentially because of the progressive 
fall in the a plot. This suggests that a consistent (i.e. reasonably scale independent) 
determination of a s is indeed going to be possible. 
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However, in the small x region, i.e. when two large scales are present (Q 2 and the 
centre of mass energy of the partonic scattering process s = ^f^Q 2 ) a NLO computation 
may be performed in a variety of inequivalent expansion schemes, which correspond to 
different ways of organizing the perturbative expansion ||. The standard two loop com- 
putation, where only logs of Q 2 are considered leading, is one such expansion scheme, 
whereas other schemes are obtained by also including logs of 1/x. Unlike different factor- 
ization or renormalization schemes, which are equivalent up to higher order perturbative 
corrections, different expansion schemes can lead to physically inequivalent predictions be- 
cause they sum up different classes of logarithms. A comparison with the data may thus 
allow one to decide in which region each expansion scheme is physically most appropriate. 

The fact that the HERA data agree well with the double scaling prediction ([!]) , which 
treats the two large scales on the same footing (retaining only terms with the same number 
of powers of logQ 2 , log K and a s ) suggests that the most appropriate expansion scheme 
in the HERA region is the double leading scheme @, where one sums all contributions to 
splitting functions where each power of a s is accompanied by either a power of log Q 2 or 
log -, i.e. the two scales are still treated symmetrically. At leading order this means that 
one should keep all terms with either one power of a s or any number of powers of a s log -, 
and at NLO all terms with an extra power of a s . The singlet splitting functions are then 
of the form@ 

P lj (x, t) = P^ (x) H |^-Pnlo( x ' *) + '") 

p&> {*) = pi 3 (*) + pi 3 (x, t) = w + pii *) , 

where (x) and P^ (x) are the usual one and two loop splitting functions, while P^ (x, t) 



1 In the standard notation in which the evolution of the singlet parton distribution functions 
f l (x,t) is given by = %*■ <8> Z 3 : different conventions were used in refs.[^j6|. 



and Pg J s (x,t) are (convergent) sums of leading and subleading singularities respectively: 

1 00 1 

P s "(x,t) = 2C A - V a n+1 -A s (t)^ n , 
x ^— ' n! 

n=3 

Pf (M) = gPf (x,t), P S 95 (M) = P**(M) = 0, 

00 1 

P«(x,t) = f ln2^T i? n / -^a n+1 -A s (t)^ 



a; * — ' n\ 

n=1 (4) 

P«/(M) = P*/0M), 



C A " ss V."") "/> 

00 



1 1 

Pf/ (x, t) = 16 In 2Ci - £ - A s (t) n £ 

n=l 
1 °° 1 

P/« (x, t) = 16 In 2C A C F - K+i -y A s (t) n C , 

n=l 



where A s (t) = 41n2— a 8 (t), a s (t) is the two-loop running coupling, and £ = log ^ for 
x < xq, zero otherwise. 

Similar expansions are taken for the singlet coefficient functions 

C*(x, t) = 8*6(1 -x) + ^ (Ci j (x) + C?(x, «))+••', (5) 

where C[ J (x) are the usual two-loop coefficient functions, while 

1 00 1 
C?{x,t) = 2 sT R n f -J2cnj^^X s (trC-\ 

n=l 

CT*(x,t) = gfCf (x,t), Cff(M) = -C™(M), Cf|(x,t) = -<7«(M). 



(6) 



The nonsinglet splitting functions and nonsinglet coefficient function remain the same as 
at two loops; there are no singular contributions in the nonsinglet channels, except in 
singular factorization schemes [JIIJ . 

The coefficients a n in (^) may be deduced PI from the BFKL kernel by means of 
appropriate factorization theorems Jl^]; a similar technique for computing a n and c n has 



been developed in 0,0,0. However, unlike the o n , the coefficients a n and c n , being 



subleading, depend on which renormalization and factorization schemes are adopted. In 



ref. [TTJJ] calculations are performed in MS, but results are also given in a corresponding 
parton scheme (DIS) in which the gluon contribution to P2 is eliminated altogether, the c n 
are zero and C i3 (x,t) = 6 %3 6(1 — x). Furthermore, the process independent singularity in 
the DIS quark anomalous dimensions, which is what makes the corresponding d n so large, 
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can be eliminated by a suitable redefinition of the normalization of the gluon distribution, 
thereby greatly reducing the coefficients a n ; this normalization (and thus this form of the 
anomalous dimension) is automatically given by off-shell factorization of collinear singu- 
larities (Qo factorization) |16[ . This factorization prescription thus defines an alternative 



parton scheme which we call QoDIS. It is also then possible to construct a QoMS scheme, 
in which the anomalous dimensions (and thus the a n ) are the same as those in the MS 
scheme, but the process independent singularity in the coefficient functions is eliminated, 
reducing the c n . Finally, it is even possible to construct a parton scheme in which both the 
a n and the c n are zero (the singular DIS or SDIS scheme) [JT7J ; all subleading singularities 
in the quark channel are then factorized into the starting distribution. 

Analytic expressions for the first fourteen coefficients a n may be found in the second 
of ref. ||12| , and for the first five a n , c n in ||14j[ , while the first fourteen a n and twelve c n in 



MS may be found in table 2 and table 3. Analytic expressions for the a n in DIS scheme 
are readily deduced from tables 2 and 3, and the corresponding expressions in the QoDIS 
scheme may be computed similarly. Numerical values for the first thirty- two a n , and a n 
in DIS are given in |J; in table 1 we give also the first twenty a n and the corresponding 
c n in MS, the c n in QoMS, and the a n in QoDIS.i 

The coefficients b n and b' n which appear in the expansion of the NLO singular splitting 
functions P§§ and P 9q have not yet been explicitly computed in any scheme. However, 
when working at NLO in the double leading scheme no computation is really necessary, 
since the requirement of momentum conservation determines these coefficients uniquely. 
Indeed, in any subtraction scheme in which there is a detailed balance of momentum 
between quarks and gluons, the splitting functions must satisfy the relations 

f dxx[P q9 (x,t)+P"{x,t)]=0, [ dxx[P qq (x,t)+P 9q (x,t)] = 0, (7) 
Jo Jo 

for all t. This equation cannot be satisfied at leading order, since P qg = P qq = 0, but at 

NLO it determined! 

bn = b' n = -a n+1 - ^fa n . (8) 



These latter coefficients are obtained [16] in practice by simply suppressing the singular 
process independent Rn factors in the expressions in ref. [13-15]. 

3 Note that if these coefficients were computed explictly in a given renormalization and fac- 
torization scheme the result would not necessarily agree with eq. (||) and thus would not conserve 
momentum at NLO in the double leading expansion. However, it is always possible to then per- 
form a change of factorization scheme such that momentum conservation is restored at NLO, and 
thus eq. (ph is satisfied. 



5 



This procedure can be extended to any finite order: for example at NNLO P£f , P 9q , P qq s 
and P qq s would be computed in MS, and then Pjjfj a and P 9q s would be fixed by momentum 
conservation. Notice that at NLO this prescription automatically preserves the colour- 
charge relation 

Pf/(x,t) = gfPf/(x,t); (9) 

it is then easy to check explicitly that because of the colour-charge relations between P 9q 
and P", P qq and P qg the supersymmetry relation 

P qq (x : t) + P 9q (x, t) = P q9 (x, t) + P"(x, t), (10) 

holds when Cp = Ca = ^Tr- 

Whereas the double leading expansion scheme discussed so far is appropriate to the 
HERA region, it is clearly inadequate at large x: there, in order to improve the two- 
loop expressions of splitting functions, it would be necessary to compute the full three 
loop contributions P3 and C 2 J ', while to include only those higher order contributions 
which have small x singularities is meaningless. Below a reference value xq the ordinary 
loop expansion scheme should thus be used instead, so the singular contributions P*- 7 , P^ s 
and C % J should all vanish there. Imposing that at x = xq the splitting functions and 
coefficient functions should be continuous, while for x < xq the logarithmic singularities 
are built up asymptotically as the evolution length xq/x increases, then leads naturally 
to the expressions (f|) and We will assume for simplicity that the parameter xq 

which sets the transition to the small x region is independent of t. We then treat it as 
a free parameter, to be eventually determined empirically. Naively we would expect it to 
lie in the range 0.01 ~ x$ ~ 0.1, at least for all reasonable schemes; here we will thus 
consider the extreme case xq = 0.1, the opposite extreme being given in effect by the 

4 Note that in [||.|6| we defined the separation of the evolution region and introduced different 
expansion schemes in different regions at the level of the evolution equation, by implicitly evolving 
the parton distribution functions at large x to give the boundary condition at x = xq from which 
they were then evolved using an expansion scheme appropriate to small x. Here the separation is 
done instead at the level of the evolution kernel, thus making the separation into large and small-x 
regions on the basis of the evolution length (in x). The present prescription has the advantage 
that it may be used in Mellin space, since the evolution remains translationally invariant. The two 
prescriptions however are identical by construction when x > xq, and only differ by subleading 
corrections in £ when x < xq. 
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two loop computation (for which xq = 0). Notice that this parameter will in general be 
factorization and renormalization scheme dependent, and should for instance be smaller 



in MS (which has large values of c n and a n ) than in the QoMS scheme. 

Since the expansions (§]) and @ are convergent, they may in practice be truncated 
after only a finite number of terms: when working in Mellin space arbitrarily accurate 
results may thus be obtained by truncating the corresponding series expansions of the 
anomalous dimensions || (even though these series have only a finite radius of convergence, 
diverging for small N). Defining 7 ZJ (iV, t) = f Q dxx N P lJ (N; t) term by term, we thus have 



7 f (N,t) ~ 2C a Xq N a n \ s {t) n - l /N n , 

7 f (N,t) = ^ 7 f (N,t), 7 f (N,t) = 7 f (N,t) = 0, 

"max 

1 q s 2(N,t) ^ fhi2C A T a n f Xo N £ a n \ a (t) n - x / N n , (11) 

n=2 

"-max 

72? (N,t) ~ -Wln2C A Xo N ^ (a n+1 + ^a^A^t)- 1 /^, 

n=2 

7L 9 (^, t) = gf 7s 9 f (iV, t), 7 f s 9 (iV, = §f 1 9 J(N, t). 

The renormalization group equation may then be integrated explicitly, while the inverse 
Mellin transform is straightforward because the anomalous dimensions ( |TT| ) are meromor- 
phic in iV: in particular there are no extra singularities or branch cuts Jj| .@ Care must be 
taken to ensure that all the subleading contributions are consistently linearized in this so- 
lution (even across thresholds). It is important to notice that this applies to the full NLO 
contribution, and not only to its large-x part; the effect of the spurious sub-subleading 
terms which would otherwise be produced may in practice be quite large. 

In order to actually perform the computation of F% we use a suitable generalization 
of the general procedure developed for two loops calculations in ref. |§. We take a set of 
parton distribution functions A = {/(x,4GeV 2 ) : / = v , q, g, . . . , x > 10~ 2 } fitted to all 
available data with Q 2 *Z 4GeV 2 , x <^ 10 -2 which predates the HERA data and is thus 
unbiased by it: in practice we work with the two MRS distributions DO' and D-' |T^ . 
We then evolve these distributions to some new scale Q 2 ,, and refit the resulting parton 
distribution functions {f(x, Qq) : / = v, q,g, . . . , x > 10~ 2 } to functions of the same form 



A quite complicated set of branch cuts are present in the analytic continuation of the complete 



series from large AA [18| . 
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as the original distributions, but with the small x behaviour of the singlet sea and the glue 
distributions now fixed to be x x for some A: 

xq(x, Ql) = A q x Xq (l - x)*fc(l + e q x 1/2 + 7„x), 

(12) 

xg{x, Ql) = A g x x ° (1 - x)^ (1 + e g x 1/2 + -y g x). 

We thus obtain a new set of distributions [A](Qq, X) which when evolved up to larger 
values of Q 2 will give (almost) precisely the same results for x ^ 10 -2 as the original 
distribution A (perturbative evolution being 'causal' both in t and l/x) but different 
results for x 10 -2 . The functions [A](Qq, A) can actually be fitted to parton distributions 
{f(x, Ql)} defined either in the MS or in a parton scheme; possible differences in the results 
will then depend on the approximation introduced by the choice of the specific functional 
form of eq. (Of) . 



Numerical computations can now be performed in Mellin space using a straightfor- 



ward extension to n max loops of the standard two loop solution [2D| of the renormalizaton 
group equations, with the form eq. (|ll]) of the small- x contributions to the anomalous 
dimensions. We have checked explicitly that the evolved parton distributions converge as 
n max is increased: very accurate results in the HERA region may be obtained using only 
the coefficients given in table 4, for all values of xq. Due to the universality of double scal- 
ing we expect the results to be insensitive to the detailed form of the starting distributions; 
the only significant constraint imposed by the large- x data is the overall normalization of 
the small-x tail. The data can thus be fitted by simply adjusting the parameters Qq, X q 
and A g , taken within suitable ranges. 

Let us first consider the two loop case (or, equivalently, double leading schemes with 
xq = 0) in the MS scheme. We begin by assuming for simplicity X q = X g = X. Still, the 
set of free parameters (Qq, X) is somewhat redundant. To understand how this works, we 
present in fig. 2a a x 2 contour plot in the A — Qq plane. All the 46 ZEUS and 76 HI 
bins which pass the cut a, p > 1 are included in the computation of the x 2 ; the cuts 
are imposed in order to only include small x data. The normalization of each experiment 
is also fitted, in accorfance witfi the expe rim ental uncertain^. Clea riy , a good fit can 
be obtained for a variety of values of the starting scale Qq; the valley floor then defines A 
as a function of Qq (or alternatively, of course, the other way around). As noted in 0, 



A useful discussion of the appropriate treatment of normalization uncertainties may be found 

in 
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reasonable results at two loops may be obtained when A(Qo) vanishes at Qo ~ 1 GeV; 
starting at a higher scale the best fit value increases, being around A = —0.2 if Qo — 2 GeV, 
consistent with the most recent MRS global fit [B21. What really matters however is not 
so much the precise value of A, but that the data prefer soft starting distributions, which 
necessarily give rise to double scaling in the HERA region: as explained in || if the starting 
distribution is too hard, say A(l GeV 2 ) ^ —0.2 or A(4 GeV 2 ) ^ —0.4, double scaling is 
spoiled and a satisfactory fit to the data is no longer obtained. Similarly the fit gets worse 
if one attempts to describe the data using a power at a much higher scale. 

In table 4 we show thus the best fit value of the parameter A obtained in the MS 
scheme at various values of Qo- Whereas the double scaling prediction of ref. with 
Qo = 1 GeV already provides a good overall description of the data, the fit improves in 
a statistically significant way when the full two loop calculation is performed, especially 
if, consistent with the x 2 plot of fig. ^|a, we choose a somewhat larger starting scale. The 
improvement in the fit is also apparent in the corresponding two loop curves in the scaling 
plot fig. [1} Evolution in DIS scheme at two loops gives almost identical results. 

We now proceed to NLO calculations in the double leading expansion scheme with 
xq = 0.1, and in the various factorization and renormalization schemes discussed above. 
The result found in the SDIS scheme (which differs from the two loop computation DIS only 
for the inclusion of gluon leading singularities 7f 9 (|TT|)) are effectively indistinguishable 
from those obtained at two loops@ (at least in the range of present-day data); we will thus 
not discuss this case further here. If the QoMS scheme is used instead the best-fit value 
of the parameter A changes appreciably; however the impact on the quality of the fit does 
not appear to be statistically significant. If we use instead the double leading MS scheme, 
the agreement with the data of the present fit is significantly worse. This, however, turns 
out to be a byproduct of the specific way we are parametrizing the boundary conditions. 

In order to understand this, we present in table 5 the results of fits analogous to 
those discussed so far, but where now we fit the small- x exponents of the quark and gluon 
distributions (|T2"D, X q and X g as two independent parameters. Furthermore, we fit the 
functional form of eq. ([1^) to parton distributions in either MS or parton schemes; in 
order to disentangle effects due to the fitting from the true scheme dependence we then 
evolve both sets of distributions in both MS and parton schemes. It then turns out that 

7 This is partly because the first few a n vanish, but also because the others are very small, so 
that their effect is concentrated in a very narrow wedge close to the boundary 
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in some schemes (and in particular in the double leading MS scheme) a good fit is only 
possible if the values of X q and X g are substantially different. In particular, at two loops 
with the fitted distribution in MS X q ~ X g , but if the distribution is fitted in DIS X q 
is significantly lower than X g . In the double leading expansion the situation is reversed: 
X q ~ X g in DIS, but X g is very significantly lower than X q in MS.! However, if the double 
leading calculation is performed with Qo factorization, then X q ~ X g in either case. Notice 
that this is genuinely an effect of the fitting, in that all of the above remains true regardless 
of whether the evolution is performed in MS or in a parton scheme. 

We conclude that in general the results of the double leading calculations do not dif- 
fer significantly from those obtained at two loops. This means that, while the leading 
singularities P/ 9 have essentially no effect at all, most of the effect of the subleading sin- 
gularities P s 9 / on F2 in the HERA region may be absorbed into the boundary condition, 
generally reducing A and increasing Qo- The HERA data on F% thus do not as yet actually 
allow us to fix the value of Xq. However, the boundary conditions themselves do change 
significantly as we go from two loops to the double leading scheme; furthermore, in the 
double leading scheme they display very strong scheme dependence. This is of course the 
essence of factorization: since the boundary condition at Q 2 , is essentially nonperturba- 
tive and scheme dependent, it must always be fitted in a comparison with the data, and 
furthermore it may change substantially when different factorization and renormalization 
schemes are adopted. This however has also the practical consequence that we may only 
reasonably take X q = X g = A in the fitted distribution if we fit MS distributions at two 
loops, and DIS distributions in the double leading scheme, unless we use Qo factorization, 
in which case we may use either. In order to avoid artifical complications related to the 
fitting procedure we shall thus henceforth perform two loop computations (fitting and, for 
definiteness, evolution) in MS and double leading computations in parton schemes; this 
will allow us to take X q = X g = X throughout. 

We can now explore in more detail the Qo dependence of A in the double leading DIS 
scheme, as summarized in the x 2 contour plot fig. ^b. The best fit A at a given scale 
is significantly reduced in comparison with the two loop case fig. ^|a because F 2 grows 

8 It is interesting to notice that the best fit X g is then rather large; nevertheless scaling is not 
spoilt. This agrees with the results of ref. ||, where we show that although at two loops only 
\ g <^ —0.3 leads to scaling, in the double leading scheme any X g <j —X s (Qo) is sufficient. 
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more quickly near the boundary as Q 2 is increased.! The scale Qo at which a given 
(soft) boundary condition should be set is thus correspondingly increased ||: the results 
obtained then differ very little from the two loop fits, except at large p i.e. very close to 
the boundary itself, as displayed by the pertinent curves in fig. [I]. Furthermore, because of 
the more significant role played by perturbative evolution, the double leading fit is more 
sensitive to the choice of Qo- However, the fit is less sensitive to the precise value of 
A. This somewhat surprising result was anticipated in ref. || where it was shown that 
actually double scaling is obtained from a wider set of initial conditions if a full double 
leading computation is performed than would be the case at two loops. 

It is now relatively straightforward to determine a s from the HERA data on F2: we 
simply include an extra parameter Aqcd (or equivalently a s (Mz)) in the fit. The results 
of such two parameter fits, for the various NLO calculations described above are displayed 
in table 6. The experimental uncertainty in the resulting value of a s (Mz) turns out to be 
remarkably small, of the order of a few per cent. The value of the starting scale Qo is seen 
to have little effect on the value of a s {Mz), though A(Qq) varies rather strongly when Qo 
is below 2 GeV (see fig. fj). 

The dependence on Q and A is explored in detail in the x 2 plots in the two planes 
perpendicular to that of fig. the a - A plane fig. ||, and the Qo-a plane fig. |j. The well- 
defined minimum in A and a s {Mz) is apparent from fig. |3]a; large values of a s are excluded 
much more strongly when the subleading singularities are included fig. ^|b. The flatness 
of the valley floor in fig. [| reflects that in fig. [2|, but again the acceptable range of Qo is 
significantly reduced when the subleading singularities are included in the calculation. In 
general all double leading calculations seem to be significantly more sensitive to the values 
of the parameters; even though they cannot yet be favoured on statistical grounds they 
could thus potentially lead to a firmer determination of a s . For the time being, however, 
we determine the central value of a s from the fits with the best x 2 , namely, the two loop 
calculations. An average of these (see table 6) gives a central value a s (Mz) = 0.120. 
The uncertainty related to the possibility of choosing various expansion schemes is then 
asymmetric, since double leading schemes tend to lead to a reduction of a s ; from table 6 
we estimate it to be +0.002 - 0.006. 

9 Very low values of the starting scale Qo (as favoured by |23| l, for example) would seem to 
be excluded in these schemes, because their effect when a s (Qo) is large is too great to admit a 
satisfactory fit to the data even with valence-like input (large positive A). 
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We now examine all other sources of error in turn. The sensitivity to the form of the 
starting distributions is explored by considering different choices of large x distributions (for 
example the two 'extremal' MRS distributions DO' and D-' [p.9| ), and also more importantly 
by varying the way in which the small-x tail is fitted to these distributions at Qq. In 
particular, choosing A s 7^ X g has little effect on a s (see table 7). Two alternative fitting 
procedures have also been tried: one in which a small-a; tail of the form x x is added by 
hand to the singlet distributions xq(x, Qq), xg(x,Qo) before the refitting (which is then 
performed over the extended range x > 10 -4 ), the normalization of the tail being fixed 
by continuity, and a similar one in which the tail has instead the form x A (l ± e-y/x), the 
parameters e q and e g being chosen such that the derivatives of the two distributions were 
continuous at the matching point (see again table 7). Estimating the error from all these 
sources rather conservatively at ±0.003, and combining it in quadrature with an error of 
±0.003 from the two parameter fit, and ±0.002 from the choice of Qq (see table 6), we feel 
confident in quoting an overall experimental uncertainty of ±0.005. 



The most important theoretical error is, just as in the large x determinations ]23 ] , the 
uncertainty due to unknown higher order perturbative corrections. This is manifest both in 



the scheme dependence (results contrasting MS and DIS schemes are presented in table 8), 
and the dependence on the choice of factorization scale M 2 = k\Q 2 and renormalization 
scale /i 2 = k 2 Q 2 (table 8). Since at present there is insufficient data to rule out all but very 
large changes in these scales empirically,0 we simply varied them by a factor of two either 
side, as is usual in the jet determinations H; this gives an uncertainty of ±0.008 — 0.004. 
It could be argued that larger variations are already ruled out in deep inelastic scattering 
by the large x data p3fl . 

The position of the quark thresholds is varied from \m q to 3m q (with m c = 1.5GeV, 
nib = 4.8GeV), to give an overall uncertainty on a s (Mz) of ±0.002: varying the prescrip- 
tion used for taking the running coupling across the thresholds gives a further uncertainty 
of ±0.001. To estimate the effect of higher twist corrections, and possible mass effects at 
the quark thresholds, we perform the analysis with a Q 2 cut on the data, discarding all 
data with either Q 2 > Q 2 ut , or Q 2 < Q 2 uV The results are displayed in table 9, and show 
a remarkable level of stability as Q cu t is varied. We can thus conclude that higher twist 
effects, since they would vary rapidly with Q 2 , must be extremely small throughout the 
HERA region. Similarly any mass effects at the quark threshold must also have only a 



10 With the exception of significant reductions in the factorization scale. 
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small effect on a s : most of the statistical significance of our determination comes from the 
data points with the smallest values of x, and thus values of Q 2 above the charm threshold 
but below the beauty threshold. 

All these various possible sources of error are summarised in table 10: combining the 
theoretical errors in quadrature, we quote a value 

a s (M z ) = 0.120 ± 0.005(exp) ± 0.009(th). (13) 

This is quite consistent with other determinations, in particular those from fixed target 



deep inelastic scattering data P^| , and the world average [Q. It is difficult to say at 
precisely which value of Q 2 the result (|13|) is obtained, since the two loop running over 



quite a wide range of Q 2 was a crucial ingredient in the match to the data @, and thus 
to the determination of a s (Q 2 ). However it should be clear from table 9 that the most 
important Q 2 range is 4GeV 2 £ Q 2 £ lOOGeV 2 . 

To summarise, we have shown that perturbative QCD works extremely well in the 
kinematic region explored so far at HERA. This remains true even when the perturbation 
series is reorganised to sum all leading and next-to-leading large logarithms, though the 
data are as yet insufficient to show to what extent this is necessary. Furthermore, since 
the shape of Fi is then relatively free of nonperturbative uncertainties, varying rather 
rapidly with a 8 , it can be used to provide a surprisingly precise determination of a s . We 
have performed such a determination using the 1993 HERA data, with the result (|T3"|). 
The experimental error, while already relatively small, should improve substantially when 
the 1994 data become available, and the expected fivefold increase in statistics might also 
provide an empirical determination of the as yet unknown parameter xq, and enable a 
reduction in the dominant theoretical error coming from the renormalization scale depen- 
dence. Although the present determination has much in common with jet determinations, 
it has the advantage that it is sufficiently inclusive that there are no uncertainties from 
hadronization, and in fact here higher twist effects seem to be very small indeed. Also, 
the summation of large logarithms is relatively easy to control, essentially due again to 
the inclusive nature of the measurement, but also to some accidental cancellations. This 
means that a further reduction in the theoretical error will presumably be possible with 
a calculation of the three loop anomalous dimensions for deep inelastic scattering, and of 
the subleading singularities at small x. 
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Figure Captions 

Fig. 1. Double scaling plots of the 1993 HERA data on Ff : i) R F F% vs. a for a) the ZEUS 
data |l|] , b) the HI data ; ii) the same vs. p. The data plotted are those included 
in the fits but with the further cut Q 2 < lOOGeV 2 ; Rf is accordingly evaluated 
with n f = 4 and (as in §) Q = lGeV, x = 0.1 and A = 260MeV. The dotted, 
dashed, dotdash, and solid curves correspond respectively to double scaling, two 
loops, double leading QoDIS, and double leading DIS, with Qq = 2GeV, xq = 0.1, 
and the remaining parameters fitted as in table 6. The curves in the cr-plots have 
p = 2.2: those in the p-plots a = 1.7. 

Fig. 2. Contours of equal x 2 in the A-Qo plane with a s (Mz) = 0.120: a) at two loops, 
MS scheme b) double leading DIS scheme with xq = 0.1. Qq is in GeV. The % 2 
is computed in the same way as in table 4. 

Fig. 3. As fig. ||, but the a - X plane, with Q = 2 GeV. 

Fig. 4. As fig. ||, but now the Qq - a plane, with a) A = —0.23 in the two loop calculation 
and b) A = —0.06 in the double- leading calculation (see table 6). 
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Table 1: Numerical values of the first twenty coefficients a n and c n in the Qo schemes and 
MS schemes. The a n in QoMS are the same as those in MS by construction, while the a n in 
DIS were given in ref. ||. 
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Table 2: Exact expressions for the first fourteen coefficients a n in MS scheme. 
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Table 3: Exact expressions for the first twelve coefficients c n in MS scheme. 
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Table 4: The fitted normalizations to the two experiments (HI and ZEUS, with uncertain- 
ties of ±5% and ±4% respectively plf), the fitted values of A and the associated x 2 s (f° r 
121 degrees of freedom). All data points with a, p > 1 are included in the fits: this cut ef- 
fectively excludes all points not inside the double scaling region. Statistical and systematic 
errors for each data point have been simply combined in quadrature. The different cases 
considered are a) the simple one loop double scaling calculation of ref. ||, with A fixed 
(Qo = 1 GeV); b) two loops (MS scheme), Qo = 1 GeV, Qo = 2 GeV and Qo = 3 GeV; 
c) two loops, DIS scheme, and double leading, SDIS scheme, both Qo = 2 GeV; d) double 
leading QoMS scheme, Qo = 2 GeV and Qo = 3 GeV; e) as d), but MS scheme. All entries 
b)-e) are full NLO computations with a s (Mz) = 0.120; the double leading calculations have 
xo = 0.1. 
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0.03 ±0.16 


64.5 




96% 


102% 


-0.26 ±0.02 


0.12 ±0.17 


72.6 




95% 


102% 


-0.23 ±0.05 


0.14 ±0.10 


57.8 


d) 


95% 


101% 


-0.25 ±0.02 


0.02 ±0.05 


62.1 




96% 


103% 


-0.26 ±0.02 


-0.07 ±0.10 


68.1 



Table 5: As table 4, but allowing X q ^ X g , and with various fitting and evolution proce- 
dures. The different cases considered are: a) MS distributions evolved in MS; b) MS distri- 
butions evolved in DIS; c) DIS ditributions evolved in DIS; d) DIS distributions evolved in 
MS. For each entry the three cases correspond respectively to two loops, double leading Qo 
schemes, double leading, standard schemes. All calculations are done with a s (Mz) = 0.120; 
the two loop calculations have Qo = 2.0, and the double leading calculations Qo = 3 GeV 
and xo = 0.1. 





norm. 


A 


a s {M z ) 


x 2 




93% 


100% 


0.42 ±0.56 


0.1183 ±0.0092 


61.4 


a) 


95% 


102% 


-0.23 ±0.03 


0.1215 ±0.0017 


58.4 




96% 


103% 


-0.26 ±0.02 


0.1212 ±0.0017 


57.8 


b) 


94% 


99% 


0.04 ±0.10 


0.1186 ±0.0043 


64.2 




97% 


103% 


-0.12 ±0.04 


0.1200 ±0.0026 


59.9 


c) 


93% 


100% 


-0.08 ±0.05 


0.1144 ±0.0034 


62.5 




91% 


98% 


-0.26 ±0.02 


0.1106 ±0.0019 


65.7 


d) 


93% 


99% 


-0.06 ±0.06 


0.1140 ±0.0040 


64.9 




91% 


98% 


-0.26 ±0.02 


0.1098 ±0.0040 


67.8 



Table 6: As table 4, but with a s (Mz) now included in the fit: a) two loops, with Qo = 
lGeV, 2GeV, and 3GeV respectively; b) double leading, QoMS scheme with Qo = 2GeV and 
3GeV; c) as b), but QqDIS scheme respectively; d) as b), but DIS scheme. 





norm. 


A 


a s (M z ) 


x 2 


a) 


95% 


102% 


-0.23 ±0.03 


0.1215 ±0.0017 


58.4 


b) 


96% 


102% 


-0.23 ±0.01 


0.1224 ±0.0018 


57.8 


c) 


96% 


103% 


-0.26 ±0.03 


0.1224 ±0.0016 


58.1 




95% 


101% 


-0.22 ± 0.03 


0.1211 ±0.0018 


59.1 


d) 


95% 


102% 


-0.23 ±0.03 


0.1212 ±0.0018 


58.7 




94% 


101% 


-0.20 ± 0.03 


0.1202 ±0.0021 


59.7 



Table 7: As table 5, but to demonstrate insensitivity to large x distributions and the fitting 
of the small x tail: all calculations at two loops with Qo = 2GeV and using a) DO', b) D-', 
c) DO', but now with X q = A, and X g = A + AA, AA = —0.3,0.3, d) using the two alternative 
fitting procedures described in the text. 





norm. 


A 


a s (M z ) 


x 2 


a) 


95% 


102% 


-0.23 ±0.03 


0.1215 ±0.0017 


58.4 


b) 


96% 


103% 


-0.25 ±0.03 


0.1215 ±0.0017 


58.2 




94% 


101% 


-0.42 ± 0.03 


0.1222 ±0.0014 


71.7 


c) 


96% 


103% 


-0.13 ±0.04 


0.1216 ±0.0022 


58.7 




98% 


104% 


-0.08 ±0.05 


0.1219 ±0.0024 


60.9 




95% 


102% 


-0.14 ±0.05 


0.1178 ±0.0024 


60.8 


d) 


95% 


102% 


-0.20 ±0.04 


0.1187 ±0.0024 


58.7 




95% 


102% 


-0.25 ±0.03 


0.1252 ±0.0018 


58.7 




95% 


102% 


-0.27 ±0.03 


0.1295 ±0.0024 


59.5 



Table 8: As table 5, but to show the scheme dependence: all calculations at two loops with 
Qo = 2GeV and using a) MS evolution, b) DIS evolution, c) as a), but with factorization 
scale ki = |,2,4, d) as a) but with renormalization scale ki = j, |,2, 4. 





norm. 


A 


a s (M z ) 


x 2 


10 


95% 


102% 


-0.24 ±0.04 


0.1213 ±0.0015 


45/92 


5 


94% 


102% 


-0.23 ±0.04 


0.1211 ±0.0024 


22/53 


4 


96% 


101% 


-0.22 ±0.06 


0.1225 ±0.0031 


11/33 


3.5 


98% 


100% 


-0.24 ±0.07 


0.1220 ±0.0043 


6/20 


3 


101% 


99% 


-0.33 ±0.13 


0.1142 ±0.0177 


3/7 


3 


95% 


103% 


-0.22 ±0.05 


0.1219 ±0.0023 


54/111 


4 


96% 


103% 


-0.21 ±0.10 


0.1221 ±0.0034 


48/85 


5 


96% 


103% 


-0.27 ±0.21 


0.1216 ±0.0039 


44/77 


7 


99% 


101% 


-0.39 ±0.35 


0.1211 ±0.0064 


26/53 



Table 9: As table 5, but to show the stability as the data are cut in Q: in the upper half of 
the table data above Q C ut are removed, while in the lower half data below Q cu t are removed. 



two-parameter fit 


±0.003 


variation of Qo 


±0.002 


fitting to large- £ p.d.f. 


±0.003 


higher order singularities 


+0.002 
-0.006 


position of thresholds 


±0.002 


form of running coupling 


±0.001 


higher twist corrections 


±0.001 


scheme dependence 


+0.008 
-0.004 



Table 10: Errors in the determination of a s (Mz)- 




Fig. 1 




Fig. 1 




Fig. 3 




Fig. 4 



